Add Butterfly Factorization #97
Merged
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
|
@Shreyas-Ekanathan can you push your version? |
| @@ -1,81 +0,0 @@ | |||
| using BenchmarkTools, Random | |||
src/butterflylu.jl
Outdated
| 🦋2!(view(Bu2, 1 : Mh, 1 : Nh), U₁u, U₁l) | ||
| 🦋2!(view(Bu2, Mh + 1: M, Nh + 1: N), U₂u, U₂l) | ||
|
|
||
| #Bu1 = spzeros(M, N) |
Project.toml
Outdated
| LoopVectorization = "bdcacae8-1622-11e9-2a5c-532679323890" | ||
| Polyester = "f517fe37-dbe3-4b94-8317-1923a5111588" | ||
| PrecompileTools = "aea7be01-6a6a-4083-8856-8a6e6704d82a" | ||
| SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf" |
There was a problem hiding this comment.
is the sparsearrays dep necessary?
src/butterflylu.jl
Outdated
| using LinearAlgebra: Diagonal, I | ||
| using LoopVectorization | ||
| using RecursiveFactorization | ||
| using SparseArrays, SparseBandedMatrices |
There was a problem hiding this comment.
sorry about that, fixed now
|
@ChrisRackauckas I think this is ready to merge. Test failures are pre-existing (someone probably should look into them though) |
|
The infos are just negative what they should be? |
|
yes, but it's not clear to me whether that's true for all of them or why they would have switched. |
|
Rebase onto latest master |
* Update README.md * Fix negative info values for NoPivot() on Julia 1.11+ Julia 1.11 changed the convention for LU factorization info field when using NoPivot(). Now, when a zero diagonal is encountered during unpivoted LU factorization, the info field should be negative to distinguish it from pivoted factorization failures. This commit updates RecursiveFactorization to match Julia's convention: - In _generic_lufact!, return negative info when Pivot=false and zero diagonal found - In reckernel!, handle negative info values correctly when adjusting offsets Fixes JuliaLinearAlgebra#95 🤖 Generated with [Claude Code](https://claude.ai/code) Co-Authored-By: Claude <[email protected]> * Update ci.yml * Update ci.yml * Add VERSION check for Julia 1.11+ negative info convention Only use negative info values for NoPivot() on Julia 1.11 and later. Earlier versions should continue using positive info values. - Added NOPIVOT_NEGATIVE_INFO constant that checks Julia version - Use this constant to conditionally apply negative info values - Ensures backward compatibility with Julia < 1.11 🤖 Generated with [Claude Code](https://claude.ai/code) Co-Authored-By: Claude <[email protected]> * Update README.md --------- Co-authored-by: ChrisRackauckas <[email protected]> Co-authored-by: Claude <[email protected]> Co-authored-by: Christopher Rackauckas <[email protected]>
ChrisRackauckas
merged commit ed51f4e
into
JuliaLinearAlgebra:master
19 of 21 checks passed
Closed
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
3 participants
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
builds on the previous PR to implement a full butterfly factorization method